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Motivated by a recent experiment by Bergeal et al, we reconsider incoherent pair tunneling in a 
cuprate junction formed from an optimally doped superconducting lead and an underdoped normal 
metallic lead. We study the impact of the pseudogap on the pair tunneling by describing fermions in 
the underdoped lead with a model self-energy that has been developed to reproduce photoemission 
data. We find that the pseudogap causes an additional temperature dependent suppression of the 
pair contribution to the tunneling current. We discuss consistency with available experimental data 
and propose future experimental directions. 



PACS numbers: 74.50.+r, 74.72.Kf, 74.40.-n 

I. INTRODUCTION 

Upon lowering the temperature, the superconducting 
gap in underdoped cuprates evolves smoothly from an 
energy gap already present in the normal stated Even 
after decades of debate, the nature of this 'pseudogap' in 
the normal metallic regime of the underdoped cuprates 
still remains a puzzle,^ and new experiments are needed 
to shed light on its nature. 

One such experiment was proposed by Janko et al.™ 
The experimental set up consists of a junction formed 
by a superconducting and a normal metallic lead in the 
pseudogap phase, separated by a tunneling barrier. If the 
pseudogap is due to the presence of preformed Cooper 
pairs, the current- volt age (/ — V) characteristics of such 
a junction should show characteristic signatures due to 
pair tunneling that differ from the standard result based 
on gaussian fluctuations.^ 

The proposed experiment was recently done by Bergeal 
et a/.p and their data appear to be consistent with gaus- 
sian fluctuations. However, even if preformed pairs is 
not a correct description, the pseudogap, regardless of 
its origin, still affects the fermions in the normal metallic 
lead, and thereby the gaussian result should not hold. A 
similar observation has recently been made in the con- 
text of the Nernst effect in the pseudogap phase of un- 
derdoped cupratesP Current vertices calculated within a 
model that rep roduces photoemission data in the pseu- 
dogap phase^ show an additional temperature depen- 
dence, which suppresses the Nernst signal relative to the 
gaussian result, consistent with experimental data. 

In the present paper, we study if a similar effect of the 
pseudogap changes the I —V characteristics in the above 
mentioned tunnel junction. We compare to experimental 
findings by Bergeal et al. and discuss possible further 
directions to improve the understanding of the incoher- 
ent pair tunneling in the pseudogap phase of cuprates. 
Throughout the paper we set h = 1 and ks = !• 



II. FLUCTUATING PAIR TUNNELING 

A direct experimental test of pairing fluctuations above 
T c is the second order Josephson effect which has 
been observed in conventional superconductors^ and 
more recently in cuprates!^ The effect is exhibited in a 
junction involving two leads, in the cuprate case, one un- 
derdoped (UD), the other optimally doped (OD), with 
critical temperatures such that T^ D < T < T° D . The 
rigid pair field of the optimally doped superconductor 
then plays the role of the external field in a typical (lin- 
ear response) susceptibility measurement. 

The net effect is that the fluctuating pairs produce an 
additional contribution 7 pa i r to the current-voltage char- 
acteristics of the junction that is directly proportional 
to the imaginary part of the pair susceptibility x of the 
pseudogap lead, 

I pah .(V,H)^eC 2 X "(q(H),oj(V)). (1) 

Here the frequency oj(V) = 2eV is linear in the bias volt- 
age V, and momentum q{H) is linear in the in-plane 
magnetic field H . The magnitude of the pair contribu- 
tion to the tunneling current is controlled by the vertex 
C which describes the pair transfer between the leads and 
depends on specifics of the junction. A measurement of 
the excess current I pa ir as a function of V and H allows 
one to trace the frequency and momentum dependencies 
of the fluctuating pair susceptibility. 

A. Microscopic theory 

In a microscopic calculation, the lowest order pair con- 
tribution to the tunneling current arises in fourth or- 
der perturbation theorypE^diagrammatically depicted in 
Fig.[l] Assuming a particle-particle i-matrix for the pseu- 
dogap lead with pairing instability in the d-wave channel, 
and keeping only the relevant d-wave part of the tunnel- 
ing matrix element, the incoherent pair contribution to 
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FIG. 1: Incoherent pair tunneling contribution, / pa i r = 
2etr (F A V AB G B G B t B G B G B V BA F A ), to the tunneling cur- 
rent. Here the trace includes summations over momenta and 
frequencies, double and single lines correspond to Gor'kov 
functions F A (ie„) of the optimally doped (OD) lead A, and 
single particle Greens functions G B (ie n ) of the underdoped 
(UD) pseudogap lead B, respectively, circles represent one- 
electron tunneling matrix elements Vj^p^and the wavy line 
denotes the particle-particle t-matrix for the pseudogap lead 
B (here for q=0). To arrive at Eq. |2| in the text, we as- 
sumed a pairing instability in the d-wave channel t B k , q (u;) = 
£q (lj) cos(2(y3k) cos(2yv) and kept only the relevant d-wave 
part oc Vi in the harmonic expansion of the tunneling ma- 
trix eleme nts (|Vpl s | 2 ) = \V Q \ 2 + \Vi\ 2 cos(2y> p ) cos(2v3 k ), see 
Refs. I5IT51 



the tunneling current is of the form (see Fig.[T]for details) 

H ) = ^Sa 2 C 2 x"{q(H),u{V)) (2) 

where x"(q, oj) = Im£q (uj) with C R the retarded compo- 
nent of fluctuating pair propagator, V and H the applied 
voltage and in-plane magnetic field, and S and a 2 the 
junction area and the lattice spacing, respectively. The 
vertex 



e„ p,k 



x cos(2(/j p ) cos 2 (2<y9k) 



(3) 



describes the tunneling of an incoherent pair^ and sets 
the magnitude of the pair current. We assume that the 
c-axis is perpendicular to the junction, and all momenta 
and coordinates contain only two-dimensional in-plane 
components. Here 



A r 



(4) 



denotes the anomalous Gor'kov function of the supercon- 
ducting lead A, with e n (fermionic) Matsubara frequen- 
cies, A p = Acos(2(^ p ) with ip p = arctan(p a /p x ) , and 
Gk(ie n ) is the single particle propagator of the pseudo- 
gap lead B (specified below). Finally, 7 = ni\Vi\ 2 /N 2 
with N the number of sites in a layer and 7ij the number 
of impurity scattering sites per unit area of the insulating 
junction. \V\ \ is defined in Fig. [ip 

The precise form of the pair susceptibility \ varies de- 
pending on the particular scenario adopted to describe 



the pseudogap phase. In this paper, we will adopt the 
standard gaussian form for the pair propagator 



1 



1 



Ao e — iau + r/q 2 



(5) 



Here, e = (T - T c )/T c , N the density of states, a = 
ir/8T and ij = ttD/8T where D is the diffusion constant. 
Alternate forms, where a is complex (in a preformed pairs 
scenario due to a BCS -BEC crossover between diffusive 
and propagating pair d 3 ! 15 ^) seems to be ruled out by 
Bergeal et al.W 

However, even if the pseudogap is not due to preformed 
pairs and a gaussian approach (Eq. 5) is relevant, the 
pseudogap will still affect the tunneling current through 
the vertex C. A similar observation has recently been 
made in the context of the Nernst effect in the pseudo- 
gap phase of underdoped cupratesP Calculation of the 
current vertices within a modeP^ which reproduces pho- 
toemission data in the pseudogap phase leads to an addi- 
tional T-dependent suppression of the Nernst effect rel- 
ative to that predicted by the gaussian model, consis- 
tent with experimental data. We will now determine if 
a similar modification occurs for the tunneling current, 
independent of whether the pseudogap is due to pairing 
or not. Before doing so, we recall that within the GGq 
approximation employed by Janko et aLpl the vertex is 



C~ —rn\Vi\ 2 N A {Q)N B (Q) 



(6) 



with Na{0) and Nb(0) the single-particle density of 
states per spin per site for superconducting lead A and 
pseudogap lead B, respectively. 



B. Pseudogap vertex 

We next investigate implications of the the pseudogap 
for the pair tunneling. Following Refs. 151171 we calculate 
the vertex ^ using the Greens function 



G k {ie n ,A B ) = 



(iei,n ~ 6c)(«eo, n + £k) - A| k 



(7) 



which is based on a phenomenological self-energy describ- 
ing photoemission data in the pseudogap phaseP^ Here 
Ak = A cos(2<pk) is the momentum dependent pseudo- 
gap, £k are the single particle energies measured from 
the Fermi level /i, and (i — 0, 1) 6i <n = e n + Ti sign(e„) 
with To the inverse pair lifetime proportional to T — T c 
(i.e., e/a), and T\ the single-particle scattering rate. For 
r = T\ = r , Eq. ([7]) reduces to the single lifetime model 



G k (ie„, As) 



U, 
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A 2 



(8) 



Eq. {8} gives a good description of the T dependence of 
the Fermi arc if T oc TP 



To compute the vertex |3]) we first derive that the mo- 
mentum sum for the two lifetime model is 

A^ 1 ^G k ( l e„)G_ k (-ie„) cos 2 (2^ k ) 

+ [f + Zl n ] K(X n ) - Y^Z hn Zl n U (Y n , X n ) } 
_N B {0) 



" A B 
where we introduced 

z . -h 



4A| 



B 



A 2 



(9) 

(10) 

(11) 
(12) 



where fc 2 - 1/(1 + e 2 /A\). 

The vertex is then obtained by completing the Mat- 
subara sum, 



37T 2 

C=—n i \V 1 \ 2 N A (0)N B (0)A(T,T ,T 1 ,A A ,A B ) (14) 



Here we have included the numerical prefactor for later 
convenience, and the temperature dependent, dimension- 
less function 



a(t, r 0s ri,AA,A B ) 

= ^M BC (T ) A A )MP s (T ) r ,r 1) A B ) (15) 



and K (z), E(z), and H(w,z) are the complete elliptic 
integrals of the first, second and third kind. 

The momentum sum over the Gor'kov Greens function, 
on the other hand, is 



iY- 1 ^^ c (ze„,A A )cos(2^ p ) 



--N A (0)k n [ E(k n ) + [E(k n ) - K{k n )] 



:N A (0)M sc (T,A A ) 



(13) 



It may then be verified that approximating in A 
the elliptic functions by their value at zero ar- 
gument describes well the temperature dependence 
of Eq. (15) (see below). Using that the ellip- 
tic functions at zero argument take the value n/2 
we may thus approximate A(T, r , Ti, A A , A B ) — 
Ai(T/Ab,To/A b ,Ti/A b ,A a /A b ) where the tempera- 
ture dependent function, normalized as Ai(0, 0, 0, 1) = 1, 

(z n = En/As) 



Al(x ,X 1 ,X2,X 3 ) 



4x 3 x 



E 



(z n + xi)[2z n + x 2 + xi] + 1 



1 /2 

{[zl + x 3 ] [(z„ + x 2 )(z n + Xl ) + ^T^ 2 + 1)] ) [(z n + x 2 )(z n + Xl ) + 1] 



(16) 



In the single lifetime model, we may 
introduce the corresponding functions 

B(T, T, A A ,A B ) = A(T, r, r, A A , A b ) and analogously 
Bi(T/A B ,r/A B , A A /A B ) for Eq. ^ in the same 
limit. As can be verified (see below), taking the 'zero-T' 
limit in order to convert the sum over Matsubara fre- 
quencies into an integral gives a good description of the 
vertex in the single lifetime model. We may thus further 
approximate B(T,T, A A ,A B ) ~ B (T/A B , A A /A B ), 
where 



B (xi,x 3 ) 



4x9 



dz ■ 



2{z + Xl ) 2 + 1 



3tt Jo ~ ^z 2 + x 3 [(z + Xi) 2 + 1] 3/2 

(17) 

Comparing then to Eq. ([6| in the previous section, 
we find that taking into account the pseudogap within 



the two and single lifetime models ^ and ^ re- 
sult in a renormalization of the pair contribution to 
the tunneling current by the T dependent functions 
.A 2 (T,r ,ri,A A ,A B ) and B 2 (T,T, A A , A B ) (their ap- 
proximations A 2 and Bq, respectively). 



So far we have neglected any external magnetic fields. 
Variation of an in-plane field allows one to trace the mo- 
mentum dependence of the fluctuating pair propagator, 
see Eq. |2|. If the relevant field scale on which the pair 
contribution gets suppressed is small enough so as to not 
affect the single particle contribution, the magnetic field 
allows one to separate these two contributions to the tun- 
neling current. We will address this matter further in the 
next section. 
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FIG. 2: T-dependence of the vertex C 2 within the sin- 
gle lifetime model for the pseudogap lead, where C oc 
B(T, r, Aa, As). Here T* is the pseudogap temperature, 
A A ~ A B = A, and r/A = V3T/T* so that the Fermi 
arcs connect at T — T* . The solid line shows the exact result 
with B given in ( 15 1 and To = Ti = F, and the dashed-dot 



line shows the approximation Bo given in ( 17 1 



FIG. 3: T-dependence of the vertex C 2 within the 
two lifetime model for the pseudogap lead, where C oc 
A(T,T ,Ti,Aa,Ab)- Here Aa — Ab = A and param- 
eters similar to Ref. 6, Ti = 200meV, A — 50meV, and 
To = (16/tt)(T - T c ), twice F GL , where T c = 80A". The solid 



line shows the exact result A given in (151 and the dashed-dot 
line shows the approximation Ai given in 1161). 



C. Comparison to experiment 

The recent tunneling experiment by Bergeal et alW 
using a YBCO/NdBCO junction with optimally doped 
(OD) NdBCO and underdoped (UD) YBCO was de- 
signed to test the above discussed theory. The corre- 
sponding critical temperatures, T~ D — 90X and T^ D ~ 
61K (the pseudogap temperature for the UD sample is 
T* ~ 250K), allow for a comparison to predictions in a 
range of temperatures extending over a considerable frac- 
tion of T c . Provided the vertex C changes only weakly 
with temperature, the data reported by Bergeal et al. are 
consistent with a standard model of gaussian fluctuations 
with susceptibility 



x"(g,w) oc 



(e + aw) 2 + (r/q 2 ) 2 



(18) 



The extracted pair contribution to the tunneling current 
is in good agreement with a Lorentzian with a width that 
was 1.6 times T GL = 8(T - T c )/n at T - T c = 6K and 
1.3 times at T — T c = 9K. This is in good accord with 
ARPES, where Tq was found to be approximately twice 
T GL W 

The vertex renormalization within the single lifetime 
model, B 2 (T,T, A A , A B ), is depicted in Fig. [2] for a T 
independent maximal value of the energy gap in the su- 
perconducting and pseudogap phase A^ ~ Ab — A with 
r scaling as T. The scaling factor used to describe the 
photoemission data is V / A = v 3T/T* , implying that the 
T dependence of B 2 {T,T, A A , A B ) is controlled by the 
pseudogap temperature T* (at which the spectral gap 
'fills up' in the antinodal region of the Brillouin zone). 
For the UD sample T* ~ 250K and the temperature 
range in which pair tunneling is detected by Bergeal et 
al. (~ 15K) is too narrow to observe any noticeable 
deviations from predictions of the simple gaussian for- 
mula Eq. ( 18 ). To test the predicted rapid suppression at 



higher temperatures coming from B 2 (T, T, A A , A B ) will 
require using a magnetic field to cleanly separate the pair 
tunneling current from the much larger normal tunneling 
current. 

Within the two lifetime model, the vertex renormal- 
ization depends on the ratios r /A B and Ti/A B . Since 
Po has a rapid T dependence near T c , one might sus- 
pect a stronger effect as compared to the single lifetime 
model. From ( 15 1 one finds, however, that for any reason- 



able value of Ti/A comparable to that found from pho- 
toemission, the dependence of A 2 (T, P , Ti, A^, Ab) on 
temperature is qualitatively similar to that in the single 
lifetime model, see Fig. [3] It becomes more pronounced 
only for small values of Ti/A <C 1. We note that a zero- 
T approximation similar to ( |17[ ) is not applicable for the 
two lifetime model. 

Finally, some remarks about the magnetic field depen- 
dence. The length scale which sets the Fraunhofer pat- 
tern (i.e. the field dependence) of the Josephson current 
is Z = X A + X B + d, where X = A tanh(VF/2A), with 
Xa.b and Wa/ B the a/b penetration depth and thickness 
of films A and B, and d the thickness of the barrier In 
the two limits of wide and thin films X(X, W) ~ A and 
X(X, W) ~ W/2, respectively. Referring then to the ex- 
perimental configuration of Bergeal et al., Wa = 200rtm, 
W B = lOOnm, d — 30nm, and the junction length 
L = 5000nm. For the optimally doped superconduct- 
ing film (A) Xa = lOOnm, and for the underdoped film 
(B) X B = 200nm if it were superconducting. How- 
ever, the latter is above T^ D and has no long range or- 
der, i.e. lead B is in the thin film limit X B ~ W B /2. 
Therefore, Xa — Xa tanh(H / A/2A J 4) — 76nm and Z = 
Xa + W B /2 + d ~ 156nm. The corresponding field scale 
H is then estimated from HLZ = 4>q where <po is the 
flux quantum, yielding H ~ 25 Gauss. Since the normal 
contribution to the current should not change much at a 
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field of 25 Gauss, while the pair contribution is strongly 
suppressed at this field, this allows one to distinguish 
the two contributions. This could be exploited in future 
experiments. 

In this context notice that, in contrast to YBCO, in 
the case of BSCCO (suggested by Janko et a'P), one has 
a stack of Josephson junctions. The effect due to the 
stack would be to create a new length scale Z' equal to 
the bilayer-bilayer separation,^ which is of the order ~ 
l.5nm. The resulting small length scale corresponds to a 
field scale H for the stack about two orders of magnitude 
larger than Z, meaning that the presence of a stack would 
not affect the I — V characteristics of the A-B junction on 
the field scale discussed above, and thus this complication 
can be ignored. 



et al. P We find that accounting for the pseudogap within 
a single lifetime model leads to a suppression of the pair 
contribution to the tunneling current relative to gaus- 
sian theory as the temperature is increased. Within the 
rather small temperature range tested in the experiment, 
however, this effect would not be noticeable. To deter- 
mine this would require differentiating the pair current 
from the much larger normal current, which would re- 
quire the application of a small in-plane magnetic field. 
We contrast this with the Nernst effect, where the nor- 
mal background is significantly smaller. Therefore, we 
suggest that such field dependent measurements be done 
in the future that would not only help identify effects due 
to the vertex, but also test the validity of Eq. 5 in the 
context of specific theories for the pseudogap phased 



III. SUMMARY 

Implications of the pseudogap on transport!^ and the 
Nernst effecP have been previously studied within a 
phenomenological model used to describe photoemission 
data. Here we study the implications of the pseudogap 
on the pair tunneling, as recently measured by Bergeal 
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